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PACS 61.25.Hq- Macromolecular and polymer solutions; polymer melts; swelling 

PACS 41.20.Cv - Electrostatics; Poisson and Laplace equations, boundary-value problems 

Abstract. - We consider a polymer brush grafted to a surface (acting as an electrode) and 
bearing a charged group at its free end. Using a second distant electrode, the brush is subject to 
a constant electric field. Based on a coarse-grained continuum model, we calculate the average 
brush height and find that the brush can stretch or compress depending on the applied field and 
charge end-group. We further look at an undulation mode of the flat polymer brush and find that 
the electrostatic energy scales linearly with the undulation wavenumber, q. Competition with 
surface tension, scaling as q 2 , tends to stabilize a lateral g-mode of the polymer brush with a 
well-defined wavelength. This wavelength depends on the brush height, surface separation, and 
several system parameters. 
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Introduction. — There are different ways to bind 
polymers to surfaces. Either by adsorption from solu- 
tion or grafting them onto the surface with a terminal 
group or having an adhering block in case of block copoly- 
mers. Such coated surfaces have many important appli- 
cations in colloidal and interfacial science. The polymer 
layer can change the hydrophobicity of the surface, pre- 
vent absorption of other molecules from solution and, in 
general, plays an important role in colloidal suspensions 
by preventing flocculation and aggregation of coated col- 
loidal particles [1,2]. 

A densely grafted polymer layer is called a polymer 
brush. The layer is grafted irreversibly on a solid sur- 
face by an end-group. Both neutral and charged poly- 
mer brushes have been studied extensively in the last 
few decades [3-11]. If there is no strong interaction be- 
tween the monomers and the surface, the brush proper- 
ties are mainly determined by the chain entropy. Neu- 
tral brushes, to a large extent, are characterized by their 
height that scales linearly with N, the polymerization in- 
dex [4-7]. Charged polymer brushes depend in addition 
on the charge density of the chain as well as the solution 
ionic strength [8-13]. 

In this Letter we aim at understanding another variant 
of polymer brushes having a terminal charge group, Ze, 
at their free end, where e is the electronic charge and Z 
the valency (see Fig. 1). The main advantage of having a 



charged end-group is that we can control the layer height 
and other properties by applying an external electric field 
and varying it continuously. This field stretches (or com- 
presses) the chains and is in direct competition with their 
elastic energy and entropy. Even in the absence of any 
external field, we expect the brush to be affected by the 
charge end-groups, because of their repulsive interactions. 
Indeed, the height profile depends on the charged group 
and an instability of the flat brush towards an undulating 
one may occur. 



Flat end-charged polymer brush. Let us briefly 
recall the equilibrium properties of a neutral grafted layer. 
The condition of highly dense layers (the brush regime) is 
£ -C R g , where R g is the chain radius of gyration, and I is 
the average distance between chains (Fig. la). The brush 
height, denned as the average distance of chain ends from 
the substrate, is denoted by h. In the 1970s, a simple free- 
energy was proposed by Alexander and de Gennes [4,5] to 
determine the layer equilibrium height. In the Alexander - 
de Gennes model, the brush height is taken to be the same 
for all chains; namely, the height distribution is step-wise. 
Later, in more refined theories, it was found that the free- 
end distribution is parabolic [6,7]. In the present work we 
remain within the step-wise distribution approximation, 
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Fig. 1: (a) Schematic illustration of a polymer brush with a 
terminal charge group. The polymer is grafted onto a flat and 
conducting surface at y = 0, while the other electrode at y — L 
provides an electrostatic potential difference V, with an aver- 
age electric field, E — —V/L. Each chain end-group carries a 
charge Ze, the grafting chain density is i~ 2 and the average 
brush height is h. (b) A continuum model of the brush used to 
calculate the electrostatic properties. The dielectric constant e 
has the same value throughout the gap between the two elec- 
trodes. The chain end-charges are bound to a two dimensional 
layer, with a charge density per unit area, a = Ze/l 2 . 



which is adequate enough as long as the wavelength of the 
predicted instability (see discussion below) is larger than 
the width of the chain-end distribution [6,7]. 

The Alexander-de Gennes expression for the brush free- 
energy is: 
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where the entropic "spring" constant is K = 9fcsT/ (Na 2 ), 
a is the Kuhn statistical length, N the polymerization 
index, fcgT the thermal energy, and vq — a 3 (l — 2\) is 
the excluded volume parameter depending on the Flory- 
Huggins parameter \- For a neutral brush, minimization 
of inrush with respect to the profile height h gives the 
well-known Alexander-de Gennes height of the brush at 
equilibrium Hq. It scales linearly with N: 
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Next, the end-charged brush is considered. As before, 
the chains are grafted onto the surface located at y = 
0, but the surface is taken as conducting and is held at 
potential ip — (see Fig. 1) [14]. At the other (free) end, 
the chains carry a charge of Ze. Without loss of generality 
we will take this charge to be positive, Ze > 0. A second 
conducting and flat surface at y = L is held at a potential 
tjj = V with respect to the surface at y = 0. Hence, the 
polymer brush is subject to a vertical average electric field 
E = —V/L. For V > 0, the field is compressing the chains, 
while for V < 0, it is stretching them. 

Although the chain elastic deformation is considered ex- 
plicitly, the electrostatic properties are calculated within 
a continuum model where the chains are coarse grained in 



the following way. The entire gap < y < L between the 
two electrodes is assumed to contain the same dielectric 
medium with dielectric constant e. The discrete end- group 
charges are replaced by a two-dimensional layer having a 
continuous surface charge density. We assume that the 
continuum limit is adequate enough and provides a good 
description of the system electrostatics. Finally, we treat 
the charged brush without taking into account the pres- 
ence of counter-ions. This assumption and the possible 
influence of counter-ions is discussed further below for a 
few relevant limits. 

The electric field has a jump at y = h, AE(h) ~ a/e, 
where a = Ze/l 2 is the layer charge density per unit area. 
For a typical grafting density of I ~ 100 nm, dielectric 
constant e ~ 10eo with £o being the vacuum permittiv- 
ity and Z ~ 5—10, the electric field jump is of order 
AE ~ 10 6 V/m = lV//im. Because the charged layer 
at y = h creates a discontinuity in the electric field E(y), 
the electrostatic problem is solved separately in two re- 
gions: the potential is marked as ?p a in the region below 
the charged layer, < y < h, and as ijjf, for the region 
above it, h < y < L. Solving the Laplace equation in the 
gap < y < L we get 
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where the coefficients b, c, and d arc determined from the 
boundary conditions: Va|y=o = 0, and ipb\y=L = V. At 
the charged layer itself y=h, the potential is continuous, 
ipa\y=h = ?pb\y=h, and the jump in its electric field is pro- 
portional to the charge density a: 

V a, , s ah V ah 
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The total free-energy can be written as the sum [15] 
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where inrush is the brush free energy, and the last two 
terms represent the electrostatic energy (in SI units). The 
volume charge density p is related to the surface one via 
the Dirac delta-function, p = aS(y — h). 

We first consider the case where the electrostatic inter- 
actions have a small effect on the thickness and we expand 
-Fbrush around its value at ho to second order in h — ho. 
The resulting total free-energy per grafting site is 
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Minimization of F with respect to the brush height h 
yields the equilibrium brush height h e \ for the charged 
case: 
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This expression is valid for low enough a 
Keh /e 2 and aV < KLh /£ 2 . 

Note that the brush height is compn 
even when the external potential gap 1 
electrodes vanishes, V — 0. In this cas 
the brush end-groups are attracted to th 
charges on the two grounded electrodes, 
tion with the closer electrode (at the orig 

In the opposite case of strong charge 
tions, the brush height can be much small( 
makes eq (6) invalid. In this case the bru 
dominated by excluded volume interactii 
eq (1)]. In the limit h <C L the brush hei 
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Fig. 2: A brush with an undulating height profile given by 
h(x) — h e \ + h q cos qx confined between two conducting and 
flat surfaces located at y = and y — L. 



Undulating charged brush. Be 

ing chains carry identical charges they i 
The electrostatic energy can be furthei 
chains compress or stretch alternatively, 
ing the distance between chain ends. In 
gate whether such an interfacial instabili 
sume that the brush height has a single i 
along one of the lateral surface dimensioi 

h(x) = h c \ + h q cos(qx) 

where h c \ is the equilibrium location of 
(7)], q = 2tt/X is the modulation q-mode. 
plitude. As was done above for the flat 
divided into two regions and the potenti 
y < h(x) and ip = ipb for h(x) < y < 
tial ip(x, y) satisfies Laplace's equation V 
following four boundary conditions: 
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Here h = (qh q sin qx, 1) / •v/T-j- (qh q sin qx) 2 is a unit vec- 
tor normal to the undulating interface given by h(x). The 
density <j c g(x) appearing in the above boundary condition 
is defined as: <7 c ff = a/ yl + (qh q singx) 2 and is related 
to the constant density a on the projected de- 
fined in the previous section for the flat layer. In our 
simplified treatment, the media dielectric constant e has 
the same value below and above the brush surface. Even 
within the uniform dielectric media assumption, we are 
able to show that the flat interface can be unstable. This 
demonstrates that the instability, due to charge-charge in- 
teractions, is different from other instabilities considered 
in Refs. [16-18], and related to heterogeneous (uncharged) 
dielectric materials placed in external electric fields. 

The potential within the polymer layer, ip a , and above 
it, ipb, are written as a power series in h q : ip a = 
E~oV>i n) and A = EZo4 n) (K) n - Nearly the 



Laplace equation is satisfied separately for each order n 
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that the zeroth order terms, ip^ ' and ipi^ , are the solution 
of the flat charged layer [eqs. (3) and (4)]. It then follows 
that for n > 

-0a™-* = (o-k^ exp(fcj/) + exp(—ky)j coskx 
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The leading contributions in h q , the layer undulation 
amplitude, can be examined by assuming that h q <C h c \ 
and expanding the electrostatic free energy up to order 
~ {h q ) 2 . We, therefore, limit ourselves to the first order in 
h q : tp = ip^ -j-ip^hq. Furthermore, we focus on the long 
wavelength limit, qh q <C 1, relevant to small amplitude 
modulations. 

For linear order in h q , only the first Fourier component 
k = q does not vanish, and we find 
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Expanding to second order in h q , both the k = mode 
and the second harmonics k — 2q do not vanish. 

The electrostatic energy difference A/ e i (per unit area) 
between the undulating (h q ^ 0) and the flat layer (h q = 
0) is given to second order in h q by 
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+ ah„ dx dz 
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cos qx 
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where L x and L z are the two lateral dimensions and h 
is h(x) from eq. (9). Straightforward algebraic manip- 
ulations give the final answer for the electrostatic energy 
difference per unit area of the brush in the long wavelength 
limit (qh q Cl): 
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The scaling of the last expression could have been guessed 
from the outset. The electrostatic energy is symmetric in 
h q — > — h q and, to lowest orders, is quadratic in h q . In 
addition, the prefactor a 2 je has dimensions of dielectric 
constant times electric field squared, and thus A/ e i must 
be linear in q. The derivation above gives us in addition 
the numerical factors and an extra dependence contain- 
ing trigonometric functions. These are symmetric with 
respect to the transformation h c \ — > L — h e \. 

Note that the externally imposed potential V does 
not appear explicitly in A/ e i. The external field simply 
stretches the brush uniformly, thereby increasing h c \. The 
interfacial instability is solely due to charge-charge inter- 
actions, as exemplified by the a 2 prefactor. The simple 
case of a thin isolated charged layer embedded in an infi- 
nite medium of uniform dielectric constant is obtained in 
the symmetric limit h e \ = ^L, and L — > oo. In this case 
one finds A/ el = — (cr 2 /2e) • qh 2 . 

Brush surface instability. — The brush instability 
mentioned above causes a deformation of the fiat layer and 
costs interfacial and elastic energy. We consider first the 
effect of surface tension and then comment on the elastic- 
ity. For a single q mode, the interfacial energy per unit 
area is A/ 7 = \"tq 2 h 2 . The total free energy difference 
Af = A/ 7 + A/ d is: 
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A/ has a minimum at a finite wavenumber q* as can be 
seen in Fig. 3(a), where A/ is plotted with a choice of 
typical parameter values. The value of q* can be obtained 
by solving the transcendental equation 
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In Fig. 3(b) we show the dependence of q* on h c \/L. The 
most stable wavenumber q* decreases monotonically as h e \ 
increases. 

In order to check whether the predicted instability can 
be seen in experiments, we estimate its order of magnitude 
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Fig. 3: (a) The brush free-energy Af from eq. (15) in dimen- 
sionless units [rescaled by {a 2 / e) 2 (h 2 q /^)\, as a function of the 
dimensionless wavenumber q [rescaled by a 2 /e~/]. Solid, dashed 
and circle-decorated curves correspond to rescaled brush height 
h e \(a 2 / e-f) = 0.5, 1 and 1.5, respectively, and rescaled film 
thickness L(a 2 /e^f) = 10. The location of the minimum in- 
creases as h c \ decreases, (b) The most stable dimensionless 
wavenumber q* as a function of rescaled height h e i. 



by taking the brush charge to be Z — 5, chain separation 
I = 20 nm, dielectric constant e = 5eo an d surface tension 
of the brush layer 7 = lmN/m. We use h e \ = je/cr 2 ~ 
10 nm and L = 10/i e i, and find q* ~ 2.8 TO 6 m _1 , resulting 
in an undulation wavelength A* = 2ir/q* ~ 2.2 fim. This 
is indeed the long wavelength limit which agrees with the 
various approximations we made. By further changing the 
system parameters: £, Ze and 7 it is possible to tune q* 
and adapt its value (in the micrometer range) in specific 
experimental setups. 

Discussion and Conclusions. — In this Letter we 
revisit the problem of polymer brushes. The new feature 
considered was to attach a charge Ze to the chain terminal 
free-end. The brush can be grafted onto an electrode (flat 
surface), while a second distant electrode is placed above 
the brush and provides a voltage gap V. The net effect is 
to have a controlled way to compress or expand the layer 
height, h. 

The charged brush creates an effective surface charge 
density that is localized at the y=h interface. These 
charges repel each other and also interact with the ex- 
ternal electric field. The equilibrium height h c \ depends 
on the competition between all electrostatic interactions 
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and the elasticity and entropy of the neutral chain, as can 
be seen from eqs. (7) or (8). It can lead to a compres- 
sion or expansion of the layer height with respect to the 
neutral brush case. But maybe the more interesting effect 
is the onset of an instability in the layer height. Em- 
ploying a linear stability analysis, the conditions leading 
to an instability of the uniform layer are analyzed, and a 
preferred wavenumbcr q* stabilized by surface tension is 
found. When elasticity of the polymer brush is included, 
it contributes a q 4 term in eq. (15) [19] in addition to the 
q 2 term originating from surface tension. The qualitative 
system behavior is similar, with a modified expression for 
the preferred wavenumber q* . 

The full system behavior in the presence of counterions 
is quite complicated and should be explored in a sepa- 
rate work, especially the limit of highly charged brushes. 
Here, we comment briefly on two extreme limits. In the 
first limit, the brush charge and the external potential are 
taken to be small enough so that the counterions are uni- 
formly distributed throughout the available volume, and 
behave like an ideal gas. It then follows that the electro- 
static potential depends quadratically on the direction y. 
In this approximation, the counterions do not contribute 
to the pressure difference AP across the brush end. The 
only source of this pressure difference is electrostatic and 
is due to the difference in \sE 2 between the two sides of 
the brush. We find AP = aV/L + a 2 /e - 3a 2 h/eL. 

In a second scenario, the brush charge is small, but the 
electrostatic energy of counterions in contact with the elec- 
trode, eV, is much larger than the thermal energy ksT. 
Here we find that all counterions migrate to one of the 
electrodes. However, since the above calculation assumes 
a fixed voltage gap, V, the electrodes will accumulate ex- 
tra charge to balance exactly the counterions. Therefore, 
the results in eqs. (6), (7) and (8) stay valid. Lastly, 
we point out that in a more physically feasible setup, the 
system may contain added salt [12, 13]. In this case, the 
electric field is screened and the brush ends do not "feel" 
the electrodes as long as the brush length is larger than 
the Debye-Hiickel screening length. 

It is worthwhile to mention some similarities between 
our charged brush and other two-dimensional systems of 
charges or dipolcs. A two-dimensional layer of electric 
dipolcs pointing in the perpendicular direction was inves- 
tigated [20] in relation with dipolar Langmuir monolayers 
at the water/air interface [20]. When the dipoles have a 
fixed out-of-plane moment but their inplane density can 
vary, a modulated phase in the inplane density can be sta- 
bilized with a preferred wavenumber q* . In addition, the 
dipolar free energy also scales linearly in q [20]. The sim- 
ilarity between the two systems can be understood in the 
following way. The charge displacement from their aver- 
age position at y = h c \ in our case is similar to an effective 
dipole whose moment points 'up' or 'down' with respect 
to this reference plane. Increase in the external E field in 
our system translates into an increase proportional to the 
average dipole strength in the dipolar system. 



More recently [21], a g-mode instability was found for 
an electric double layer where a charge density bound to 
a surface was allowed to fluctuate laterally. The model is 
motivated by an experimental setup where charged am- 
phiphilcs coat heterogeneously a mica surface. In the ex- 
periment the surface contains patches of positive and neg- 
ative charges but the overall surface charge (summed over 
all patches) remains zero. As the surface was placed in 
contact with a salt solution, a local electric double layer 
is formed. Positively and negatively charged counterions 
are attracted, respectively, to negative and positive charge 
domains, resembling our system as well as the undulating 
dipolar one [20,21]. 

We hope that the simple considerations mentioned here 
will motivate experimental studies of end-charged brushes 
where some of our predictions can be tested. 
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